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SUMMARY 


In  order  to  investigate  the  dynamic  stability  of  a  parachute  the 
system  of  the  non-linear  equations  of  motion  is  derived.  In  contrast 
to  previous  contributions  this  system  of  differential  equations  is 
not  linearized,  but  is  integrated  by  means  of  electronic  computers  for 
a  number  of  examples  of  personnel  guide  surface  parachutes.  The 
influence  of  the  suspension  line  length,  the  effective  porosity  of  the 
canopy,  the  weight  of  the  load  and  the  apparent  mass  on  the  dynamic 
bevahiour  are  studied  for  che  examples  chosen.  The  results  show  that 
the  oscillations  of  parachutes  cannot  be  described  by  linearized 
equations  of  motion.  For  instance,  only  the  non-linear  procedure 
enables  the  complete  description  of  the  velocity  oscillations.  Further¬ 
more,  it  becomes  evident  that  frequency  and  damping  depend  on  amplitude. 
For  some  amplitudes  the  square  of  the  oscillation  period  is  proportional 
to  the  length  of  the  suspension  lines.  With  increasing  porosity  of 
the  canopy  the  damping  increases  and  the  oscillation  period  decreases. 


SOMMAIRE) 


Afin  d’  enqu@ter  sur  la  stability  dynamique  d'un  parachute,  on  a 
^tabli  le  syst&me  des  Equations  non-lindaires  du  mouvement. 
Contrairement  A  des  contributions  pr^cfCdentes,  ce  syst&me  d’ Equations 
differentielles  n’est  par  lin^aris^,  mais  il  est  int£gr£  au  moyen  de 
calculatrices  ^lectroniques  pour  plusieurs  exemples  de  parachutes  de 
surface  guides  de  personnel.  L’  influence  de  la  longueur  de  la  ligne 
de  suspension,  la  porosite  effective  de  la  calotte,  le  poids  de  la 
charge  et  la  masse  apparente  sur  le  comportement  dynamique  sont 
^tudids  comme  exemples  choisis.  Les  re'sultats  montrent  que  les 
oscillations  des  parachutes  ne  peuvent  etre  d^crites  par  des  Equations 
de  mouvement  lin^aris^es.  Par  exemple,  seule  la  procedure  non-lin^aire 
permet  la  description  complete  des  oscillations  de  vitesse.  De  plus, 
il  devient  Evident  que  la  frequence  et  1’  amortissement  dependent  de 
1’ amplitude.  Pour  certaines  amplitudes,  le  carr^  de  la  p^riode 
d’  oscillation  est  proportionnel  &  la  longueur  des  lignes  de  suspension. 
Avec  une  porosity  accrue  de  la  calotte,  1’  amortisseiu-nt  augmente  et  la 
p^riode  d’  oscillation  d^croit. 
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NOTATION 


velocity  vector,  resultant  velocity,  velocity  components  (m/sec) 

rate  of  descent  (m/sec) 

angular  velocity  (sec-1) 

suspended  load  (=  m;g)  (kp) 

mass  of  the  load  (kp,  sec2/m) 

apparent  air  mass  (kp.sec2/m) 

ratio  of  apparent  mass  to  mass  of  the  load  (=  nig/mj)  (1) 
angle  between  parachute  axis  and  the  vertical  (1) 
angle  between  velocity  vector  and  the  vertical  (1) 
angle  of  attack  (1) 

moment  of  inertia  of  canopy  (kp.m.sec2) 
radius  of  inertia  (m) 

distance  from  load  point  to  centre  of  gravity  of  canopy 

internal  force  (kp) 

external  force  (kp) 

drag  (kp) 

lift  (kp) 

moment  (m. kp) 

coefficients  of  drag,  lift  and  moment  of  canopy  (1) 

projected  area  of  inflated  canopy  (=  rV)  (m2) 

radius  of  inflated  canopy  (m) 

damping  (sec-1) 

oscillation  period  (sec) 

coordinates  fixed  in  earth  (m) 


v 


p 


air  density  (kp.  sec2/®4) 

acceleration  due  to  gravity  (=  9.81)  (m/sec2) 


g 

Time  derivatives  are  designated  by  a  dot. 

Subscripts 

c  canopy 

l  load 

x  ,  y  in  direction  of  x-,  y-axis  of  parachute 
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INVESTIGATIONS  ON  THE  DYNAMIC  STABILITY  OF 
PERSONNEL  GLIDE  SURFACE  PARACHUTES 

R.  Ludwig*  and  W.  Heins* 


1.  INTRODUCTION  -  PROBLEMS  AND  PURPOSE 

OF  THE  INVESTIGATIONS 

The  dynamic  stability  of  parachutes  is  of  great  importance  in  many  respects. 
Insufficient  stability  of  a  man-carrying  parachute  may  lead  to  injuries  while 
landing;  a  lack  of  stability  in  deceleration  parachutes  can  strongly  disturb  the 
control  of  the  aeroplane  during  the  landing  run;  and  instability  of  a  stabilizing 
parachute  of  a  missile  or  a  space  vehicle  would  be  contrary  to  the  real  purpose  of 
such  a  parachute. 

Advances  in  aeronautical  engineering  and  the  growing  field  of  application  of 
parachutes  have  in  the  last  decades  led  to  the  development  of  parachutes  of  high 
dynamic  stability.  In  particular,  favourable  canopy  shapes  with  an  appropriate 
porosity  have  been  developed  experimentally.  Ribbon  parachutes  and  guide  surface 
parachutes  are  typical  examples  having  good  stability  characteristics. 

The  theoretical  investigation  of  the  dynamic  stability  had  been  undertaken  quite 
early;  that  is,  already  by  1918,  as  can  be  seen  in  a  paper  by  Brodetzky1.  But  in 
more  recent  times  there  have  only  been  a  few  papers  on  this  subject,  for  example  by 
Henn2,  and  more  recently  by  Riabokin3,  of  the  Department  of  Aeronautical  Engineering 
of  the  University  of  Minnesota,  Minneapolis  (USA);  by  J.F.  Reagan  and  F.J.  Stimler" 
(USA),  and,  very  recently,  by  W.  G. S.  Lester7  (Great  Britain).  All  these  papers  are 
confined  to  the  linearized  theory  and  the  assumption  of  small  disturbances.  Sample 
calculations  can  only  be  found  in  the  paper  by  Henn2. 

It  is  the  purpose  of  the  investigations  described  in  this  Report  to  calculate  a 
greater  number  of  cases  for  the  personnel  guide  surface  parachute  by  means  of 
electronic  computers  after  having  derived  the  general  non-linear  equations  of  motion. 
The  results  will  be  discussed  with  a  view  to  understanding  the  different  modes  of 
behaviour  in  dependence  on  certain  parameters,  which  will  again  be  useful  with  regard 
to  experimental  investigations.  The  prototype  of  the  guide  surface  parachute,  on 
which  is  based  the  personnel  type,  is  an  invention  during  World  War  II  of  Prof. 

Dr.  H. G.  Heinrich,  at  that  time  member  of  the  Flugtechnisches  Institut  Stuttgart 
(Germany),  under  the  direction  of  Prof.  Dr.  Madelung.  Three-component  measurements 
carried  out  by  H.G.  Heinrich  and  E.  L.  Haak  have  greatly  contributed  to  these 
investigations5. 

The  principal  reasons  which  have  led  to  these  theoretical  investigations  were: 

(a)  The  conditions  for  a  linearized  theory  -  the  existence  of  small  disturbances  - 
are  not  satisfied  in  practice. 


*  Deutsche  Forschungsanstalt  ftlr  Luft-  und  Rumfahrt  e.V.,  Institut  fSr  Flugmechanic , 
Braunschweig  33,  Flughafen,  Germany 
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(b)  The  linearization  of  the  equations  of  motion  leads  to  a  decomposition  of  the 
system  of  differential  equations;  the  equation  governing  the  velocity 
disturbance  splits  off  and  yields  a  monotone  but  not  a  periodic  decrease  with 
the  frequency  of  the  parachute  oscillations,  as  would  be  expected. 

(c)  The  curves  for  the  aerodynamic  coefficients  CD  ,  CL  ,  and  CM  in  dependence 
on  the  angle  of  attack  show  that  a  linearization  does  not  meet  the  true 
conditions,  for  3cM/Bac  sometimes  even  takes  a  negative  value  in  the 
environment  of  ac  =  0  . 

Hie  oscillation  behaviour  of  a  given  type  of  parachute  depends  on  numerous 
parameters;  the  most  important,  which  are  to  be  varied  here,  are  the  length  of  the 
suspension  lines,  the  weight  of  the  suspended  load,  the  assumptions  on  the  apparent 
mass  of  the  canopy,  and  the  porosity  of  the  parachute  cloth.  The  oscillation 
behaviour  is  characterized  by  frequency  (i.e.  oscillation  period)  and  damping,  as 
generally  known  of  linear  systems.  Detailed  calculations  show  that  discrepancies 
occur  due  to  non-linear  behaviour  and  also  that  different  initial  conditions  lead  to 
a  different  oscillation  behaviour.  All  these  results  can  be  considered  to  be  useful 
for  experimental  investigations. 

Even  though  comparative  measurements  on  the  personnel  guide  surface  parachute  were 
not  possible,  observations  on  other  types  show  qualitative  agreement  with  the 
calculations. 


2.  EQUATIONS  OF  MOTION 

To  begin  with,  the  equations  of  motion  on  which  the  calculations  are  based  will 
be  indicated.  The  detailed  procedure  is  given  in  Appendix  1.  The  equations  are 
valid  under  the  following  conditions: 

(a)  The  system  canopy  -  load  is  assumed  to  be  rigid. 

(b)  The  mass  of  air  associated  with  the  inflated  canopy  must  be  considered  an 
inertial  mass,  but  not  a  heavy  mass4.  It  is  designated  as  apparent  mass. 

The  weight  of  the  canopy  can  be  neglected. 

(c)  The  aerodynamic  forces  which  act  on  the  load  are  assumed  to  be  negligibly 
small.  This  assumption  is  not  satisfied  with  loads  of  special  shape,  as, 
for  instance,  with  recovery  systems  of  rockets  or  with  stabilization  systems 
of  torpedos. 

Figure  1  shows  the  geometric  quantities  and  the  forces  acting.  The  differential 
equations  of  forces  result  of  the  equilibrium  of  the  mass  forces  with  the  internal 
and  external  forces  acting  in  the  direction  of  the  tangent  and  of  the  perpendicular 
to  the  trajectory. 

For  the  canopy  we  have 


mc(7c 


p  (i) 
rc 


+  F„ 


(e) 


(1) 


3 


and  for  the  load 


The  equation  for  the  moments  of  the  canopy  in  the  system  of  coordinates  fixed  in 

the  axis  of  the  parachute  and  referred  to  the  load  point  is 

-  mcs\y  -  racsvcxw  =  ‘Me  •  O) 

If  we  introduce  the  condition  of  rigidity,  the  geometrical  relationship  between 
the  angles  a  ,  /?  ,  y  ,  the  external  forces  with  the  aerodynamic  forces  in  the  usual 

form,  where  the  coefficients  depend  on  the  angle  of  attack  ac  of  the  canopy,  we 

obtain  a  system  of  four  differential  equations  of  the  1st  order  (see  Pig.  2): 

(m;  +  mc)(vx  -  wvy)  +  mcsw2 


=  mjg  cosy  SVg[Cp(ac)  cosac  +  CL(ac)  sinaj  (11)* 

(m;  +  mc)(vy  +  covx)  -  mcsa> 

P  2r 

=  -  ni;g  sinry  +  —  SvclCD(ac)  sinac  -  CL(ac)  cosacJ  (12)* 

C* 

2  2'  •  ^  2 
mc(i^  +  s2)o>  -  mcs(vy  +  a>vx)  =  -  —  S2Rv*CM(a(;)  (13)* 

7  =  o).  (14)* 

The  parachute  trajectories  come  from  a  further  integration  of  the  velocity 
I  components  in  the  system  of  coordinates  fixed  in  the  ground.  Thus  the  motion  of  load 

and  canopy  is  completely  described  by  the  following  15  quantities: 

for  the  load: 


v  ,  vY  ,  vv,y,o),/3,  a,  x'  ,  y'  ; 
x  y 

and,  for  the  canopy,  as  vx  =  vcx  : 

vc  -  vcy  -  ac  -  4  •  xc  -  • 

3.  DATA  FOR  CALCULATIONS  ON  THE  PERSONNEL 
GUIDE  SURFACE  PARACHUTE 

The  calculations  on  the  personnel  guide  surface  parachut  (see  Pig. 3)  are  based  on 
the  following  dimensions: 


Numbered  as  in  Appendix 
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R  =  3. 0  m,  s  -  9. 1  m,  i*  =  0.344m2. 

To  determine  the  apparent  air  mass  the  volume  of  the  canopy  is  first  substituted 
by  an  ellipsoid  of  revolution  with  a  volume  equal  to  that  of  the  canopy.  The  centre 
C  (Fig. 1)  of  the  ellipsoid  coincides  with  that  of  the  canopy.  The  ellipsoid  is  then 
enlarged  to  take  into  account  the  apparent  air  mass.  The  moment  of  inertia  Ic  or 
the  radius  of  inertia  ic  can  then  be  determined  from  known  relations.  For  a  usual 
type  of  man-carrying  parachute  with  a  normal  load  (%  100  kp)  it  can  be  assumed, 
according  to  von  Karman6  that  the  apparent  mass  is  of  approximately  the  same  order  as 
the  mass  of  the  load,  i.  e.  u  =  m  /m,  ^  1  . 

As  initial  conditions  we  assume  that  the  parachute  which  is  steadily  descending  at 
a  rate  v0  is  deflected  at  an  angle  yc  ,  as  in  the  case  of  a  gust.  Hence,  for 

t  =  0  :  v  =  v0,  y  -  yQ,  w  =  £  =  0, 
where  vQ  results  from  the  equilibrium  condition  Wj  =  D  : 

vQ  =  ^2W,/PSCD0  . 

The  value  0.25  has  generally  been  chosen  for  yQ  ,  corresponding  to  14.3°.  The 
wind-tunnel  measurements5  which  have  already  been  mentioned  comprise  the  aerodynamic 
coefficients  Cx  ,  Cy  and  CM  referred  to  the  axes  x  ,  y  in  the  system  of 
coordinates  fixed  in  the  parachute  (see  Fig. 4).  The  coefficients  CD  ,  CL  and  CM 
(for  different  lengths  of  the  suspension  lines)  were  obtained  from  Cx  ,  Cy  ,  CM  . 

All  aerodynamic  coefficients  apply  to  a  parachute  of  the  same  geometric  shape,  but 
are  given  for  different  effective  porosities.  The  effective  porosity  V  ,  introduced 
by  H. G.  Heinrich8  is  a  non-dimensional  quantity  providing  the  ratio  of  the  mean  flow 
velocity  through  the  parachute  fabric  to  the  velocity  of  the  free  air  stream.  The 
curves  plotted  for  CM  show  at  once  that  the  parachute  is  statically  stable  when  the 
greatest  porosity  rj  =  0.096  is  used,  as  ^Cn/^ac  >  0  •  But  when  77  =  0.042  or 
77  =  0  (non-porous  canopy)  the  parachute  is  statically  unstable,  for  in  the  environment 
of  a  =  0  :  3c„/3a.  <  0  . 

C  M  C 

4.  RESULTS  AND  DISCUSSION 

Considering  first  the  statically  stable  case  (77  =  0.096),  oscillations  of  a  certain 
frequency  are  expected  to  be  damped  with  time  (Fig.  5).  The  following  is  of  interest 
here: 

The  velocities  v  ,  vc  ,  and  the  component  vx  show  a  two-fold  frequency  compared 
to  the  velocity  components  v„  ,  v„„  ,  to  the  angular  velocity  co  and  to  the  angles 

/DO  ** 

y  .  P  ,  Pc  ,  ct  ,  and  aQ  .  From  the  physical  point  of  view  this  is  easy  to  understand 
when  the  motion  of  the  parachute  is  considered  to  be  that  of  a  pendulum.  From  the 
mathematical  point  of  view  this  can  be  explained  from  the  terms  of  o>2  and  o>vy 
of  the  differential  equations,  which  must  immediately  lead  to  an  oscillation  of  a 
doubled  frequency. 
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The  snail  amplitudes  of  oscillation  of  the  canopy  compared  to  the  oscillation  of 
the  load  are  remarkable.  This  can  be  seen  from  the  angles  ft  and  /?  shown  in 
Figure  5.  This  means  that  the  motion  of  the  load  is  similar  to  that  of  a  pendulum 
with  a  length  s  . 

The  two  frequencies  of  the  velocity  oscillations,  as  illustrated  in  Figure  5,  are 
confirmed  by  cine-theodolite  measurements  of  a  parachute  trajectory.  The  measurements 
were  evaluated  in  the  D.F.L.9.  Such  a  case,  though  for  a  different  type  of  parachute, 
is  shown  in  Figure  6,  giving  the  vertical  and  the  horizontal  velocity  components  v^ 
and  v^  .  As  a  comparison  the  corresponding  calculated  quantities  are  plotted  for 
the  personnel  guide  surface  parachute. 

At  first  sight,  the  quantities  plotted  versus  time  do  not  show  a  remarkable 
difference  with  respect  to  the  oscillation  behaviour  of  a  linear  system.  They  become 
more  distinct,  however,  by  representation  in  the  phase  plane  (Fig.  7),  where  u>  =  y 
is  plotted  versus  y  as  is  generally  done  for  a  technical  oscillation  process.  The 
time  scale  of  the  curve  shows  that  the  oscillation  period  increases  with  decreasing 
amplitude.  Of  the  amplitude  ratios  it  becomes  evident  that  damping  decreases  with 
amplitude.  The  representation  in  the  phase  plane  shows  the  oscillation  period  T 
and  the  damping  k  referred  to  the  mean  time  t  of  a  complete  oscillation  (Fig. 8). 

At  the  same  time,  the  load  is  varied,  while  the  apparent  mass  is  kept  constant.  It 
can  be  stated  that  the  oscillation  period  hardly  depends  on  the  mass  of  the  load. 

The  oscillation  period  of  a  mathematical  pendulum,  for  instance,  is  independent  of 
the  pendulum  mass,  but  this  is  not  valid  any  more  for  the  physical  pendulum.  The 
variation  of  period  with  amplitude  can,  however,  be  approximately  10  per  cent.  This 
is,  all  the  same,  enough  to  be  taken  into  account  in  experimental  investigation,  at 
least,  when  appreciating  the  scattering  of  results  for  a  drop  test  disturbed  by  wind 
and  gusts.  Damping  is  rather  strongly  influenced  by  amplitude  (up  to  about  50  per 
cent),  a  bigger  load  has,  of  course,  a  stronger  damping  effect. 

The  trajectory  curves  of  canopy  and  load  (Fig. 9)  will  be  considered  next.  In  these 
curves  the  parachute  is  plotted  for  time  intervals  of  1  second.  On  the  one  hand,  the 
porosity  is  varied  (tranjectories  a,  b,  c).  We  see  that  with  increasing  porisity  there 
is  an  increase  in  damping  and  a  decrease  in  the  oscillation  period.  On  the  other 
hand,  the  length  of  the  suspension  lines  is  varied,  while  the  porosity  is  kept  constant, 
and  it  appears  that  the  oscillation  period  increases  with  increasing  length  of  the 
suspension  lines  (CDo  as  well  decreases  with  increasing  porosity). 

The  time  history  for  the  velocities  and  angles  (Fig.  10)  shows  that  a  more  precise 
statement  can  be  made  on  the  oscillation  period  in  dependence  on  the  pendulum  length; 
this,  of  course,  pre-supposes  that  the  oscillation  period  refers  to  the  same 
amplitude,  hence  to  the  initial  value  y0  .  It  can  be  seen  (Fig. 11,  upper  part) 
that  the  square  of  the  oscillation  period  is  proportional  to  the  distance  canopy-load. 
This  suggests  basing  an  empiric  formula  on  that  for  the  mathematical  pendulum,  where 
a  factor  of  proportionality  must  be  introduced.  The  calculation  of  a  number  of 
cases  for  three  different  porosities  shows  that  the  drag  coefficient  CDo  can  be 
introduced  as  a  factor  of  proportionality  and,  substituting  this  factor  by  the  rate 
of  descent,  we  obtain 


T  %  c 


Do 
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The  oscillation  periods  determined  from  this  approximate  formula  (dashed  lines  in 
Fig.  11)  agree  quite  well  with  those  of  exact  calculations.  The  damping  is  scarcely 
affected  by  the  length  of  the  suspension  lines  (Fig.  11,  lower  part);  in  the  case 
7]  =  0.096  a  flat  maximum  appears  for  the  usual  length  of  the  suspension  lines.  It 
is  open  to  question,  however,  as  to  how  far  this  statement  can  be  generalized, 
especially  when  3cM/3ac  takes  on  a  negative  value,  and  the  oscillation  becomes 
indifferent. 

Up  to  now,  the  apparent  mass  had  been  set  up  so  that  pc  =  1  for  W(  =  100  kp  . 

As  an  assumption,  the  influence  of  a  change  of  pc  (Fig.  12)  will  still  be 
demonstrated.  According  to  the  conception  gained  that  motion  of  the  load, 
essentially,  is  similar  to  that  of  a  pendulum  of  length  s  ,  it  can  easily  be 
realized  that  the  oscillation  period  is  practically  independent  of  the  apparent  mass, 
apart  from  the  already  known  dependence  of  amplitudes.  It  is  remarkable,  as  can  also 
be  seen  from  the  trajectories,  that  at  equal  initial  position  the  canopy  will  be  more 
deflected  with  the  smaller  apparent  mass  on  account  of  its  smaller  inertial  forces; 
at  larger  apparent  mass,  of  course,  the  damping  will  become  larger.  The  amplitudes  of 
load,  on  the  contrary,  will  be  enlarged  with  increasing  pc  only  in  a  very  small 
degree. 

In  all  cases  considered  up  to  now,  the  air  density  p  =  pQ  on  the  ground  was  the 
basis  of  the  calculations  and  it  was  kept  constant  for  the  calculated  section  of 
trajectory.  Now,  cases  with  different  air  densities  or,  in  other  words,  with 
different  altitudes,  will  be  considered,  but  the  respective  air  density  within  the 
calculated  section  of  trajectory  is  again  assumed  constant. 

The  following  points,  however,  are  pre-supposed: 

(a)  the  aerodynamic  coefficients  remain  constant,  i.e.  the  porosity  does  not  change; 

(b)  since  the  apparent  mass  decreases  at  lower  air  densities,  we  assume  that  the 
apparent  air  volume  is  kept  constant,  which  means  that 


The  subscript  o  relates  the  respective  values  to  the  altitude  0  . 

Comparing  the  time  histories  (Fig.  13),  say  /3C  ,  for  air  densities  corresponding 
to  the  altitudes  between  0  and  10  km  ,  it  can  be  seen  that,  in  the  first  instance, 
the  same  phenomena  occur  as  when  varying  the  apparent  mass.  Since  also  the  aero¬ 
dynamic  forces  decrease  in  correspondence  with  the  altitude,  the  effect  will  be 
enlarged.  This  means  that  the  lateral  deflection  is  larger  than  for  the  same  p.  on 
the  ground;  also,  the  damping  is  somewhat  larger,  and  the  oscillating  period  smaller, 
with  increasing  altitude. 

In  a  linear  oscillating  system,  the  character  of  the  oscillation  is  not  influenced 
by  initial  conditions,  whereas  in  a  non-linear  system  this  problem  is  much  more 
complicated.  This  is  particularly  significant  in  those  cases  in  which  no  restoring 
moment  occurs  in  the  environment  of  ac  =  0  (Fig.  14).  Starting  from  a  small  initial 


deflection,  this  deflection  increases  at  first  until  the  canopy,  on  account  of  the 
oscillating  load,  inverses  its  direction  of  movement.  Figure  14a  shows  that  this  is 
the  reason  for  the  lateral  drift  of  the  parachute,  Starting  from  a  fairly  large 
deflection  (Fig. 14b),  an  almost  indifferent  oscillation  without  lateral  drift  is 
obtained  on  account  of  the  restoring  moment. 


5.  CONCLUSION 

Of  the  numerous  examples  which  have  been  treated,  it  becomes  evident  that  parachute 
oscillations  show  many  typical  modes  of  non-linear  oscillation,  thus  confirming  that 
the  true  behaviour  of  a  parachute  is  not  correctly  represented  by  a  linear  procedure. 
Phenomena  such  as  the  dependency  of  amplitude  on  frequency  and  damping  cannot  be 
ignored  in  connection  with  experimental  investigations.  By  experiment  on  other  types 
of  parachutes  the  occurrence  of  two  frequencies  of  the  ratio  1:2  could  be  qualitatively 
confirmed  as  well  as  the  order  of  magnitude  of  these  frequencies.  It  had  already 
been  supposed  that  the  relation  between  oscillation  period  and  length  of  suspension 
lines  could  well  be  understood  by  application  of  the  pendulum  formula,  but  up  to  now 
this  could  not  be  derived  from  exact  calculations.  The  indicated  empiric  formula  can 
also  be  obtained  under  certain  simplifying  assumptions  from  the  system  of  the 
linearized  equations.  The  results  of  the  calculations  show  that  the  effect  of  the 
apparent  mass  is  not  very  important. 

As  it  principally  affects  damping,  it  should  be  rather  difficult  to  obtain  more 
accurate  estimations  from  drop  tests  which  are  never  carried  out  in  undisturbed  air. 
That  the  initial  conditions  are  of  special  interest  for  the  time  history  of  the 
trajectory  curves  is  particularly  easy  to  understand  with  regard  to  the  complicated 
relation  between  the  moment  coefficient  and  the  angle  of  attack.  In  practice  it 
seems  to  be  of  importance  that  under  certain  initial  conditions  a  drift  may  occur 
for  the  same  parachute,  while  this  would  not  be  the  case  under  different  initial 
conditions. 

The  authors  feel  these  investigations  may  make  some  practical  contribution  to  the 
problem.  Whether  the  assumptions  made  are  satisfied  in  practice  should  be  examined 
by  test. 
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APPENDIX  1 
Equations  of  Motion 

We  introduce  the  internal  forces  Pc^  and  F(^  which  act  between  canopy  and 
load,  and  the  vector  of  the  angular  velocity  <o  (perpendicular  to  the  vertical  plane 
in  which  the  motion  takes  place). 

The  equations  of  forces  are  considered  separately  for  the  centre  C  of  the  canopy 
volume  and  the  load  point  L  .  Therefore  the  tensions  in  the  suspension  lines  F  ^ 

•4  /  j  \  C 

and  Fj  '  act  as  internal  forces.  The  following  vector  equations  yield  the 
equilibrium  of  the  mass  forces  with  the  internal  and  external  forces  in  the  direction 
of  the  tangent  and  of  the  perpendicular  to  the  trajectory: 

=Fc(i)+Fc(e)  (A.  1) 

m;(v+  [3v])  =  .  (A. 2) 

To  these  equations  we  add  the  moment  equation  referred  to  the  load  point,  viz. , 

-  ">cs(ycy  +  vcxw)  =  "Mc  •  (A-3> 

the  vectors  being  split  in  components  of  the  system  of  coordinates  fixed  in  the  axis 
of  the  parachute. 

As  the  system  canopy-load  is  assumed  to  be  rigid,  the  relation 

Fc(i)  =  P;(i)  (A. 4) 


must  be  applied. 

For  the  velocities  we  can  deduce: 


v 


cx 


V 


X 


vcy  '  vy  =  ’sw  • 
Figure  1  gives  the  angular  relations 

7  -  P  +  a  =  Pc  +  ac  ] 

tan  a  -  -  -L 
vx 

tan  a  =  -  .  cy  . 

vv 


(A.  5) 
(A.  6) 


> 


(A.  7) 
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x'  =  v  vos/3  (A.  15) 

y'  =  -v  sin/3  .  (A.  16) 

Furthermore,  to  carry  out  the  calculation,  nine  algebraic  relations  are  necessary: 


’ey 


Vy  -  SW 


4 


2  2 
V  +  V 
x  T  vcy 


(A. 17) 

(A.  18) 


= 


arctan  f — 


a  = 


-  arctan/ — 
V. 


4  =  r  - 


P  =  y  -  a 

Xp  =  x'  -  s  cosy 
y£  =  y'  +  s  sin  y 


(A.  19) 

(A.  20) 

(A. 21) 

(A.  22) 

(A. 23) 
(A.  24) 
(A.  25) 


Equations  (A. 11)  to  (A.  25)  are  used  for  the  complete  description  of  the  behaviour 
of  canopy  and  load;  by  integration  of  the  system  under  given  initial  conditions  one 
obtains  in  dependence  on  time: 


for  the  load:  v  •  vx  •  vy  >  "X  >  a>  ,  ft  ,  a  ,  x'  ,  y' 

and,  for  the  canopy:  vc  ,  vcy  .  4  •  ac  -  xc  •  yc  • 

If  we  substitute  cos  aQ  =  vx/vc  and  sin  ac  =  -vcy/vc  in  Equations  (A. 11)  and 
(A. 12),  and  if  these  two  equations  are  solved  for  vx  and  vy  ,  while  vy  in 
Equation  (A. 13)  is  replaced  by  the  corresponding  expression,  we  obtain: 


vx  =  o)Vy  -  CjW2  +  C2  cosy  -  C3vc(vxCd  -  vcyCL)  ,  (A. 26) 

vy  =  wvx  -  [C,  siny  +  C5vc(vcyCD  +  vxCL)  +  C6v2CM]  ,  (A. 27) 

"  =  -  [C7  siW  +  CeVvcyCD  +  vxCL)  +  C9VcCM]  •  28> 


with  the  coefficients: 


12 


C 


1 


Mcs 


pS 


2m  ( 

(1 

+  MC) 

PS(12 

+  s2) 

2m 

<N 

sg 

N 

p 

1  + 

pSR 

pc  m,N  ' 


=  “g(i  *  +  s2) 


pSRs 

m,N 


pSs 
2m  jN  ’ 


(A.  29) 


N  =  (1  +pc)i2  +  s2  . 


Equations  (A. 26)  to  (A.  28),  together  with  Equations  (A. 14)  to  (A. 25)  now  represent 
the  system  on  which  are  based  the  sample  calculations.  Hie  integration  of  this 
system  of  differential  equations  has  been  carried  out  on  a  digital  computer  as  well 
as  on  an  analogue  computer.  The  numerical  operations  have  been  programmed  according 
to  the  integration  method  provided  by  Runge-Kutta  for  the  digital  computer  Siemens 
2002  in  the  Centre  of  Calculation  of  the  D.  F.L.  At  steps  of  0.1  second  for  an 
integration  range  of  sometimes  60  seconds,  the  required  accuracy  was  reached.  Only 
every  second  step  was  printed.  Parallel  operations  were  carried  out  on  an  analogue 
computer  PACE  231  R  in  the  D.  F.  L. 
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fjr' 


Fig.  1  Symbols  and  equations  of  motion  I 

(mL  +mc I  fv„-u)VyJ  +mc  sco2 

=mtg  cos  y  -  jSv2[Cd  (uc)  cos  ac + CL  (a,.) sin  aj 
(ml+mc)(Vy+u)vx)-mcsu> 

=mtg  sin  y  +  j-S v2[C0 (ac)sin  ae-CL( ac) cos  aj 
me(ic2+s2)w-mcs(vy+(o  vj 

=-js2Rv2CM(ac) 

ysiCD 

Fig.  2  Equations  of  motion  II 


Fig.  3  Personnel  guide  surface  parachute 


see  HG  Heinrich  and  £  L  Hook,  ASD-TDR’62-100  (!962) 


Fig.  4  Aerodynamic  coefficients 
0:  t)  -  0;  □:  rj  -  0.  042;  A:  ri  =  0.  0 


Pig.  5  Time  histories 
W*  =  100  kp;  y0  =  0.25;  17  =  0.096 


Calculated  results 
(personnel  guide  surface 
parachute) 

Wt  =  100  kp 
yo=0.25 
1/  =  0.096 


Experimental  results 
obtained  from  cinethe- 
dolite  measurements 
(flat  circular  parachute) 
Wt  =  !00kp 


based  on  a  paper  by  0.  Weber { 
WGL- Meeting  1962 

20  22  24  26  26  30sec32 

t 


Fig. 6  Comparison  between  calculated  and  experimental  results 
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Pig. 11  Oscillation  period  T  and  damping  k 
W j  =  100  kp 


Pig. 12  Angle  between  the  velocity  vector  and  the  vertical  of  the  canopy,  the 

apparent  mass  being  varied 

Wj  =  100  kp  ,  70  =  0.25  ,  V  =  0.096 


20 


- H=  0  km ,  /ic=  WOO 

-  2  km,  0.622 

-  6  km,  0.538 

.  10  km,  0.338 

Pig. 13  Angle  between  the  velocity  vector  and  the  vertical  of  the  canopy,  the 

altitude  being  varied 

W,  =  100  kp  ,  7„  =  0.25  ,  7?  r  0.096 


Fig.  14  Trajectories 

Wj  =  100  kp  ,  17  =  0  ;  (a)  J0  =  0. 10  .  (b)  70  =  0.40 


DISTRIBUTION 


Copies  of  AGARD  publications  may  be  obtained  in  the 
various  countries  at  the  addresses  given  below. 

On  peut  se  procurer  des  exemplaires  des  publications 
de  1* AGARD  aux  adresses  suivantes. 


BELGIUM  Centre  National  d’  Etudes  et  de  Recherches 

BELGIQUE  A^ronautiques 

11,  rue  d’Egmont,  Bruxelles 


CANADA 


DENMARK 

DANEMARK 


FRANCE 


GERMANY 

ALLEMAGNE 


GREECE 

GRECE 


ICELAND 

ISLANDE 


ITALY 

ITALIE 


LUXEMBURG 

LUXEMBOURG 

NETHERLANDS 
PAYS  BAS 


Director  of  Scientific  Information  Service 
Defense  Research  Board 
Department  of  National  Defense 
'A'  Building,  Ottawa,  Ontario 

Military  Research  Board 
Defense  Staff 
Kastellet,  Copenhagen  0 

O.N.E.R.A.  (Direction) 

25,  Avenue  de  la  Division  Leclerc 
Chdtillon-sous-Bagneux  (Seine) 

Zentralstelle  fur  Luftfahrt- 
dokumentation  und  -information 
Miinchen  27,  Maria-Theresia  Str.  21 
Attn:  Dr.  H.J.  Rautenberg 

Greek  National  Defense  General  Staff 

B.  MEO 

Athens 

Director  of  Aviation 
c/o  Flugrad 
Reykjavik 

Ufficio  del  Generale  Ispettore 
del  Genio  Aeronautico 
Ministero  Difesa  Aeronautica 
Roma 

Obtainable  through  Belgium 


Netherlands  Delegation  to  AGARD 
Michiel  de  Ruyterweg  10 
Delft 


NORWAY 

NORVEGE 


PORTUGAL 


TURKEY 

TURQUIE 


UNITED  KINGDOM 
ROYAUME  UNI 


UNITED  STATES 
ETATS  UNIS 


Mr.  0.  Blichner 

Norwegian  Defence  Research  Establishment 
Kjeller  per  Lillestrom 

Col.  J. A.  de  Almeida  Viama 
(Delegado  Nacional  do  'AGARD' ) 

Direc$ao  do  Servi90  de  Material  da  F.  A. 

Rua  da  Escola  Politecnica,  42 
Lisboa 

Ministry  of  National  Defence 
Ankara 

Attn.  AGARD  National  Delegate 

Ministry  of  Aviation 
T.  I.  L. ,  Room  009A 
First  Avenue  House 
High  Hoi  born 
London  W. C. 1 

National  Aeronautics  and  Space  Administration 
(NASA) 

1520  H  Street,  N.W. 

Washington  25,  D. C. 


$ 

Printed  by  Technical  Editing  and  Reproduction  Ltd 
Harford  House,  7-9  Charlotte  St.  London.  If.  1. 


CO 

CO 

CO 

CO 

a 

o> 

CO 

CD 

CO 

T3 

T3 

O' 

lO 

O' 

in 

«H 

<H 

CO 

CO 

CO 

CO 

o 

o 

CO 

CO* 

CO 

CO 

CO 

CO 

lO 

IT) 

s 

to  Q 
A  U  t? 

„  o  o 

<5  £ 


O. 

•*h  c_> 

»H  * - •  . 


+>  Jl  o  <M 

OJ  O 

S  a  5  £ 

•  0  fc  O 

4)  .  o  >rm* 

s8s  « 


P  l  ;  .  ©  0) 

o  a  «  ©  .a  ja 

a  o  fl  "o  p  p 

a  a  o  s 

3  h  *j  o  ►, 

o  S 

j2  «  rH  O 

»  p  ^  Q  H  q 

01  »H  B  fit!  3 

p  *3  a  §  8 

43  c  g  g  *H  ^ 

s §  0  °t  I 

£  .  ►»  v  o  +* 

3  m  p  a  a  *. 

q  g  *H  ©  O 

Q  g  o  a  p  <=> 

VH  §  O  a 

n  r>  H  'O  •*> 

°  -s  41  to  P 

§ s ► s* s 

§  -H  0)  -H  *D  O  . 
^  .  a  ao  m  n 
p  £  p  e  r*  o  © 


LJ  a  T3  a  p 

8  g-S  §  8  3  S 

o  o  p  *h  3  ’o 

°  «H  ft  >,  +3  g  _ 

U  -P  -H  O  cj  £  T3 

"  O  P  G  rH  PI  O 
S  B  o  0)  .H  g  -H 

^  w  a  •*■»  .2  p 

J,  »W  0)  O'  O  ^  0) 

1  °^  £  §2  Q 

5  »S  g 

k  o  4>  +3  a  ^  -H 

S  s  ^5  £  {j 

5  s  e  **  o  .„• 


03  CD 

X3  X3 


.-h  «a 

X2  CD 

as 

s  s 


2.  o  p  0)  ©  1-1 

d  a  h  a  o 


(0 

p 

C7  w 
"  0) 

a 

© 

& 

rH 

T3 

a 

to 

© 

u  P 

p  s 

•M 

► 

© 

o* 

to 

© 

© 

M  »H 

(0 

xs 

© 

p 

e 

S  a 

s 

0 

(0 

MH  8  o 

ft)  a)  H  tO 
J3  n  ft  a> 
+J  43  |  X3 

p  p  x  a> 
0  0  #  x» 


-  ©  41 

03  x  a 

X3  P  p 

Sis 


<PJ 

0 

nc3 

©  83 

G  bo  0) 

c 

•r  1 

0 

•H 

P3 

§ 

C  p  0 
0  (0  to 
ft  d  V 

Ui 

P  ©  _ 

•H 

CJ 

c3  'D 

p 

•H  5  0 

0 

03 

w 

=3 

CL> 

C3 

0  ^  05 

T3 

03 

S5  u 

<13 

4—i 

©  g  0 

03 

a  ^  -H 

ocJ 

p  p 

C3. 

^=: 

s 

Q 

=  oi  3 

G 

©  G 

G 

03 

rH 

CO 

©  *K  P 
X  tfl 

0  <f  © 

o3 

X2 

03 

p  c  x) 

•H  p 

o> 

etS 

©  5 

H  ©  P 

GS 

G 

C=J 

»  Q  0 

£ 
a  p 

e 

•H 
r— * 

03 

03 

(x 

C-3 

03 

©  « 

■0  a  B 

3  g  to 

P  w  M 

roc© 
c  c  to 

xa 

3 

•o 

•P  4j  03 

•-*  n  0> 

41  0 

•o 

03 

G  g  P 

an  £ 

a? 

O 

03 

ro  8j  0 

X2 

O 

f— . 

3  4.  C 

3  0 

n— 4 

P 

CZ3 

«  H  » 
© 

F~> 

O 

A 

E 

41  O 

9  x  M 

flj  4)  rH 

CQ 

P< 

03 

0  c 

J3  J3  G 

a> 

t3 

03 

03 

G 

• 

to  p  ■X 

M  ft 

P  <H  X 

<L> 

*r“ < 

0  ©  a 

0  0  © 

xs 

» — < 

fc-  iH  ”0 

g  £  §  5  5  8  | 

£  —  TJ  **  3  2  © 

o  *o  ©  -  a  p 

45  ©  O  0)  b  p  o 

9  5  P  fa  ®  «h  ° 

.bO  T?  *H 

W  P  a  B  ©  ° 
oio  ms  bo 

f-i  w  u  a)  o  a 

t©  oj  a  w  ti  *h 

•OOP  Ml  G 

-  o  P  P  g  B 

x  ©  *h  p  o  ©  a 

6)  H  £  t,  H  T) 


U  P  I  . 

3  g  §  S 

.s°s 

2  0  0>  P 

3  Q)  *H 

U  p  k.  -H 

1  3  rH  O 

J  a  a  g 

2  8  °  0 
’  £  .  t>> 

i  S SS 
1  S  Ss. 
!»S> 

■I  §  — 1  OJ 

; 3 

>  “  £  p  ■ 

>  2  ° 

>  O  •  § 

*  w)  C 

!  0  *h  a 

3  ©  ’[! 

>  s  o  P 

>5  BO 
^  to 

>  -  <H  © 

;  a  o  •o 
5  ©  ©  ‘ 

1  **  a  p 

’  *  o  ©  ^ 

!  W  a  § 

S  " 

a  -a  S  ’ 

W  S  ! 

!  g  5  to  i 

'  ©  g  2  I 

esli 


<0 


P  O  P  |  .  .  O  O 

s  s § Ssa s  s 

^  n  d  H  n 

9)  «  g  H  ft  rH  2 

JT3  «  H  B  d  9 
pop^pacg 

P  g  J?  O  O  ° 

&11§§g5| 

O  W  U  >J  ’O  o  ^ 

9a)  &  Jj*  p  8  o  ^ 
op  2oak<o 
d  m  9  o  o  ft  ^ 

o  o  .2  rt  D  ►> 
A  h  £  V  M  P 

p  g  M  *2  ►  bC  *H  *H 

o  c  »  a  to 
P-Ho^^^^O 

o  >  ^  ,  fl  ft  O  l. 

Cfl  p  P  B  *H  o 
►»  43  gj  P  fc,  at*  a 
^  o  ■D  l1 

•H  Z  IJ  °  ft  tt 
®  !J  o  T3  a 

P  °  o  J  °  9  S  -h 

k«  *<“<  ro  CtH  raOy) 

O  S  o  °  +>  «H  fl 

£.  a  u  *h  a  ►,  +3  s 

p  ax  -p  *h  o  ed  J" 
a)KS°^o»Ho 

►  •OSBOOrHS 
P  ♦*  CO  3  •H  » 

P  O  ,,  P  O  0*0*^ 
oxasoTJOtn^ 

o  p  2  h  0  5 

P  5  M  0)  *W  ^ 

(H  C  ^  Q  P  HI  *H 

oj  .  o  *H  'rv  5  ** 

s  3  s  1 1 5  s  „■ 

-  B  m  5  o  P  0)  O 
43  So  fl  h  a 

PpP  OOWP 
bo  C  "  T3  XJ  T5  3  *H 

So  2  ^  p  p  o' 
h  ”  8  ►  w  g 

h  a  *>  2j  «  a)  o 

ft  **  T  0)  o  p 
o  ft  5  «h  to  43  « 
B(d2g£4»*Jfl 
•h  iP  ?  d  G  a) 
hoi^RoqmX 
43  i  0>  O  0)  w 

fl  p  o  o  2 

°*o  c  K  J  J  J  w 

a  g  g**gsa  • 

©•S-S^gEl^ 

go  *H  j_,  O 

0)  O  h  B  _ 

OCDr-HWkdOlO^ 

xi  x  no  d  a  to  p 
PPETSOP.  Off 
.  a  fi  H  k<  fl  E 

p  p  k  o  *h  p  o  o  a 

OOOJjQrHfcfcrHT 


)  p  I  .  .  O  O 

J  o  fl  w  O  43  J3 

>  a  o  g  *o  p  p 

fl  fl  O  3 

!  a  ■h  p  o  ►, 

O  O  P  p  p  s 

<  a  d  h  o 

i  S  *3  a  §  S 

fl  ®  SI  ci 

I  a  °  t  5 

i  £  ,  >,  o  o 


a  0 

0  a 

ft  y 
p  § 

P 

O 

O 

O  O 
ft  ki 
O  ft 

0 

0 

43  k* 

O  P 
«  w 
S  0 
O  ■n 

H 

O 

*o 

to 

>> 

p 

P  3 

0 

P  fH 

► 

O 

tfl  *H 

5^ 

•H 

to 

0 

0  5 

■H  1 

43 

ft  O 

k. 

to 

ofl 

p 

P 

St 

0 

0 

*  *2 

fl  P 

p 

ft 

to 

p  0 

Q 

•o  0 

(4 

;  i  §  °  §  §12  s 

i  m  +>  -h  5  a  E  -o 

f  5  I  b  S  3  "  5 

■  3  ■H  5  ►< 

^  ih  a  o'  o  ^  <u 

1  E 8 3  n 

S  m  »  <W  +J  a 

k  5  *  P  fl  *  -H 

■  ,H  rH  a  P  P 


.ISS5S  a 
!  »  g  §  *  o  a  3 

i  5.  o  p  o  o  p 

So  fl  h  fl  o 
>  o  o  3  p  to 

!  „  «n  fl  13  fl  P  O 
'3  m  P  *H  C 
'  S  S  ►  W  fl  O 

I  S  *2  w  ®  o  3 

t»  ^  ta  O  O  p  p 

f • S3  SS  S« 

e  S  g  g  „  |  g 

£0005? 

”  O  'O  2  tfl 

*  bK  9  *°  3  W  ® 

S  &  3  +>  S  3 

»  _  T5  'H  in  2  0) 

O^o  .  Q,  La 

X3  O  o  O  B  P  O 

0  5  P  u  "  im  B 

*  V  b  o  ^  -h 

to  ^  ft  B  O  0 
OO  C  B  U> 

i-h  to  h  o  o  -fl  fl 

to  od  .a  to  P  iH 

•OOP  be  fl, 
„  fl  k.  k.  fl  B 

K  O  P  fl  O  O  a 

a)  fl  H  £,  t,  H  t! 


